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INTRODUCTION 


pr«d -oj-tiveiy 

recent studv (ref rvf ^ applied to combined environments* A 

and jLdea^tSrj^^ "»°- 

coflmunity ap«cifically needs methods (n%* comfort. The ride qeelity design 

dn™t.'pi::r -- r ;^fuc“ie"?iy‘-:^s'J?y :c."oJ"th'r- 

lndl,ido.l nolU aid vibi«tirSi;Lr“«riri 

evaluating eomparative ride comtorti and ( 4 ) (or dev^loif™"*"'’ *’’ accurately 
criteria of passenger eccepteii.! levelcping accurate and verifiable 


extenirv.”’rn:«ch‘°prlL.”^°S a''”*"" •■> 

for estimating passenqir ride f"<> comprehensive empirical model 

nor noise. ?hrS’ndriL?^ri a*? * fcnctlon of vehicle vibration and inte- 
which obtained subjective comfort r^^*^ ”»* developed in the research program, 

exposed to controlled combinations of inse'i!ij"i”r«lon’irtr”°"° 
quality apparatus at NASA Lanqlev Researrvh r 4 »r, 4 .xm passenger ride- 

unique capability of transformina {rw^'iw'i /4 i resultant model has the 

environment into subjective discomfort unTt of the noise and vibration 

to produce a single discomfort index typifvLa"L^p" combining the subjective units 
ment. This capability is onrof ^ passenger acceptance of the environ- 

Results of individual studies in the approach, 

references 2 through 11, anfthe reLuf^t 

reader is referred to thesTr^^ences for ZT\ . P>^fented in reference 12. The 

for definition of the discomfort index and for^df*^ information on the model studies, 
literature. acomrort index, and for discussion of the ride-quality 


To apply the model presented in 

comfort requires spe Ific knowledae of fhf 2 f assess or predict vehicle ride 
noise and vibration. ?or examole in of the interior 

ter frequency, and bandwidths of do " ®®veral of the vibration axes, the level, cen- 

not alwuys u?kightfc^«rrd! This lx 

spectra has indicated that mode^comfort eatlLw"’ the model to actual vehicle ride 


This paper presents and describes a 
model that removes the necessity that the 
a part ot model input. The computational 


simplified version of the NASA ride comfort 
user identify vibration spectral peaks as 
procedures required to obtain discomfort 
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estimates are discussed, and a user-oriented ride comfort computer program is 
described. The simplified model is applied to sample helicopter and automobile ride 
environments . 


DESCRIPTION OP SIMPLIFIED NASA RIDE COMPORT MODEL 
Development of Simplified Model 

Because of the complexity and judgment required in identifying spectral details 
for input to the original NASA ride comfort model, an alternative method was con- 
sidered. This method Involved the concept of applying frequency weightings to the 
vibration acceleration spectra within each axis of vibration. The frequency weight- 
ing for each axis reflected human comfort sensitivity to vibration frequency for that 
particular axis. The root-mean-square (rms) acceleration level of each weighted 
spectrum was then determined and input to the comfort model. 

The frequency weighting factors for vertical and lateral vibration were based on 
human comfort sensitivity to narrow-band vibration (bandwidth of 2 Hz). The weight- 
ing factors for roll vibration were based on human comfort sensitivity to sinusoidal 
vibration. Longitudinal and pitch vibrations were weighted uniformly over the fre- 
quency range from approximately 1 to 5 Hz because human comfort sensitivity to vibra- 
tion in tl.ese axes had been obtained only for random vibration centered at 3 Hz with 
a bandwidth of 5 Hz, The weighting factors for the five vibration axes are listed in 
table I. When applying these weightings to actual vibration acceleration spectra, 
interpolation may be necessary between frequency values. The weighting factors for 
each aixis are plotted in figure 1 . 

The above frequency weighting of vibration can be justified by considering the 
vertical and lateral vibration conditions. The procedure involved computation of the 
discomfort produced by a wide range of random vertical and lateral vibrations having 
different center frequencies and bandwidths. These computations were performed using 
the single-axis discomfort eqv itions for vertical and lateral random vibration given 
in reference 12. For ease of computation, the vibrations were assumed to have the 
spectral characteristics given in figure 2. For each center-frequency and bandwidth 
combination given in reference 12 the discomfort was computed for root-mean-square 
(rms) acceleration levels of 0.02g and O.IOg. (The acceleration due to gravity, g, 
is 32.2 ft/sec^.) Next the square of the vertical or the lateral weighting function 
(figs. 1(a) and 1(b)) was applied to each spectrum (i.e,, center-frequency and band- 
width combination), and the weighted rms acceleration level was computed for each 
axis. The resultant weighted rms vertical and lateral acceleration levels were plot- 
ted against the corresponding discomfort predictions, as shown in figures 3 and 4. 

The correlation coefficient between computed discomfort and the weighted vertical rms 
acceleration levels was 0.970 with a standard error of estimate of 0,28. For lateral 
vibration the correlation coefficient was 0,948 with a standard error of estimate of 
0,19. These high correlations indicated that in the range of random vibration char- 
acteristics for which the model is valid, weighted rms acceleration is a good dis- 
comfort estimator. The final step was to relate the weighted rms acceleration to 
subjective discomfort units by computing the least squares linear regression fit tg 
the data in figures 3 and 4, The resulting equations could then be used to compete 
subjective discomfort when weighted rms vertical and lateral accelerations were 
known . 


A similar approach could not be applied to verify use of frequency weic.iting for 
roll vibrations. However, based on the results obtained for vertical and literal 
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vibrations, we assumed that the weighting function derived from the sinusoidal roll 
discomfort equations of refejrence 12 would produce similar results. 


Model Inputs 

Vibration." The basic inputs to the model are the weighted rms vibration levels 
in one or more of the five axes of vibration. The user must apply the frequency 
weighting (table I) to determine the weighted rms acceleration levels of each con- 
tributing axis. For purposes of this model, the vibration can be measured at any 
floor location of the user's choice. 

Noise.- The noise inputs required by the ride comfort model are the A-weighted 
noise levels measured at approximate head level above the point at which vibration 
measurements are made (or at a location that the user considers representative of the 
vehicle interior noise environment). The A-weighted levels used are those within 
each of six octave bands having center frequencies of 63, 125, 250, 500, 1000, and 
2000 Hz. 

Trip duration .- An optional model input is trip duration (in minutes) defined as 
the elapsed time from the start of a vehicle ride. Duration is included if the user 
wishes to correct discomfort estimates for adaptation effects. 


DISCOMFORT COMPUTATION 
Procedure 

The procedure required to compute estimates of passenger discomfort in a given 
ride environment is as follows: (1) compute discomfort due to vibration in each 

applicable axis; (2) compute combined-axes discomfort; (3) correct for trip duration 
effect, if desired; (4) compute discomfort due to interior noise in each octave band; 
(5) compute total noise discomfort; and (6) compute total estimated discomfort. The 
step-by-step computation procedure is described below. 

Step 1 . Compute discomfort due to vibration in each applicable axis .- The equa- 
tions relating discomfort to weighted rms acceleration for random vibration in each 
axis are 


vert 


'0.241 + 44.672 (g ) 

’w vert 

,68.772 (g ) 
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w lat 
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2.406 0, 
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•ft 

For 6^ > 0.10 rad/aec^ 
For 0^ < 0.10 rad/aec^.- 
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f 0.41 + 5.07 6 , 
^pi tch ■' 


.8.62 


For > 0.116 rad/sec^' 

For < 0.116 rad/sec 
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(5) 


w)iiGr6 0 Do 

DISC^ uni^rw Dpi tch are estimated discomfort values in 

<r ! Z T ' t' ' '-P»«iv,ly; 

levels ^9w)lonq» and are the weighted rma acceleration 

a unJL 1 - longitudinal, roll, and pitch axes, respectively? in 

g units for linear accelerations and in rad/sec-^ for angular acceleration. 

When the vibration is known to sinusoidal in the vertical, lateral, and roll 
axes, then the estimated discomfort is predicted from 




(Diat)sine - 


^Droll^sine = 



K, (g ) 

1 vert 

For 

'«3>v.Ct ' 

^vert 

For 

'5,’v.rt < 

■4 >9,1 

For 

'«,’lat » ‘>•0691 

i^lat 

For 

*’s’lat < 

•^7 ®s 

For 

5, > 


For 

*s < >'9l 


( 6 ) 


(7) 


( 8 ) 


where the subscript sine indicates that the estimated values of discomfort are due to 

(nr™,) 

In^ >in«,old,l vibrations in the vertical, lateral, and roll a«ee, 

SoSf ^ V t nnd rad/seo^ tor angalar aoceUra- 

tions, and Kq, Ki, ..., Kg are empirical constants given in table II for sinusoidal 
vert, .cal vibrations, in table III for sinusoidal lateral vibrations, and in table iv 
for sinusoidal roll vibrations. Longitudinal and pitch vibra?ionr;rrassSmSl ^ 51 
random in nature and are always computed using equations (3) and (5). 

the U 3 er^mu 3 t^fi??t^co?Sh'^^^h^^a^ ' ^^^^^^ compute combined-axes discomfort, 

the user must first compute the discomfort produced by combined vertical, lateral. 


.,n^^ <^iscomfort values are in subjective discomfort units called DISC. These 

ratio relationship to one another. For example, if d , « 2 and D - i ?Ln 
vibratf™! L Z Z&I ' “ 


4Tfc. <• 
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and roll vibrations and then the additional discomfort due to combined longitudinal 
and pitch vibrations. These computations are performed as follows: 

1. Rank the discomfort values for vertical, lateral, and roll vibrations that 
were computed in step 1. Let D, , D 2 , and D 3 represent the highest, middle, and 

lowest discomfort values , respectively, 

2m Compute 


», - o,/o^ 


<91 


■>4 - Fa" * 


( 10 ) 


>C1 -R 


2 2 2 
+D2 ^03 


( 11 ) 


-0.44 + 1 .65 D 


cl 


combi 


1.14 D 


cl 


For > 0.881 

For D , < 0.88 1 

Cl 


( 12 ) 


3. The discomfort, in DISC units, due to vibration in the vertical, lateral, and 
roll axes is 


4 combi 


- D^)/0.40 


D 


VI.R 



The value of DyLR 
use. 


combi 

obtained from 


equation ( 


For D. < 0,40 

4 

R^ > 3.0 
Otherwise 

13) should be retained 



(13) 


for subsequent 


4. Set D 5 equal to the larger and Dg equal to the smaller of the discomfort 
values calculated from equations (3) and (5) for longitudinal and pitch vibrations. 

5. Compute 


"2 - 


( 14 ) 


“02 ■ Ts" ♦ “e" 


(15) 


D 


comb 2 


10.70 D 


c 2 


«c2 

For 

%2 

> 

1 . 0 ^ 


For 

°C2 

< 

1 .of 


(16) 
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6, The discomfort in DISC units, due to vibration in the longitudinal and pitch 
axes is 


Dr ^ D.(D . . - D_)/0.40 
5 6 comb 2 5 


D 


LP 


D 


comb 2 


For < 0*40 and^ 

> 3.0 
Otherwise 


i 


(17) 


The value of Dj^p from equation (17) must be retained for subsequent use. 
7. Compute the combined-ajtes discomfort in DISC units, from 


D 


comb ® rVLR * °LP 


(18) 


Step 3. Correct for trip duration effect .- The NASA ride comfort studies indi- 
cated that passengers tended to adapt to vibration environments presented for dura- 
tions up to 2 hours* If desired, a duration correction, in DISC units, can be 
computed as follows: 


/t).0031 - 0.012 t 


AD. *(-0.72 
dur 


For 1 min < t < 60 min 
For 60 min < t < 120 min) 
For t < 1 min 


(19) 


where t is the elapsed trip time in minutes. Duration corrections were not deter- 
mined for trips longer than 1 20 minutes . We anticipate that fatigue may become a 
factor for very long trips and might tend to counteract the adaptation effect 
described above. No data are currently available, however, to define specific 
effects of longer durations. 

The total duration-corrected vibration discomfort is given by 


D = D . + AD . 
VIB comb dur 


( 20 ) 


Step 4. Compute discomfort due to interior noise in each octave band .- Noise 
discomfort contributions are computed for noises within six octave bands having cen- 
ter frequencies of 63, 125, 250, 500, 1000, and 2000 Hz, The model does not apply to 
noise contained within octave bands outside this range. The A-weighted sound pres- 
sure levels within each of the six octave bands are required. The model Is 

applicable for = 65 to 100 dB. If the noise level within an individual octave 
band is leas than 65 dB(A) , then that octave band is assumed to contribute zero 
discomfort. The .\ctual level of discomfort produced by any of the octave bands 
depends on the level of vibration simultaneously present in the environment. The 
equation used to estimate the noise discomfort contributed by an octave band of noise 
is 


^N(i, j) 


(A. + 
3 




)W. 

1 


(21 ) 

/ 
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where noise discomfort contribution due to a noise level of j dB(A) 

within the octave band with center Irt-quency at i Hz, in tlie presence of vibration 
that produces a discomfort level of ^nd Hj are empirically determined 

coefficients (see ref. 11) tor eacli noise level in the ranqe frr^m 65 to 100 dB(A) 
(see table V); and is a weiqhting factor (see ref. 11) that corrects for the 

effect of the octave tend at i Hz (see bible VI). Rquation (21) is applied for 
each applicable octave band. 

Step 5. Compute total nois e di seomfort If noise is present at a sufficient 
level in more than one octave band, then the total ncise discomfort is computed from 
the following equation: 

Noise ^N(i,j)max ^ ^N(i,j) ^N(i,j)max 

where DN(i,j)max discomfort associated with the octave band that produces 

maximum discomfort and implies summation over octave bands. 

Step 6. Compute total estimated discomfort The total discomfort produced by 
the combined noise and vrbrationehvironraent is given by 




^Noise ^ 


D. 


VIB 


(23) 


Discussion of Model Outputs 

The various discomfort parameters computed in the preceding section provide the 
ride development specialist with a variety of useful information for use in assess- 
ing and diagnosing vehicle ride-quality problems. The total discomfort index, 
equation (23), corresponds to the overall average discomfort level that a given ride 
environment would produce. It incorporates the combined effects of the noise and 
vibration present in the environment and can be used for direct comparison of ride 
comfort in different ride environments. When used in the design stage, it provides 
estimates of vehicle ride comfort in the presence of anticipated noise and vibration 
levels. 

One of the more powerful aspects of the NASA ride comfort model is the capabil- 
ity to separate and identify the relative contributions of noise and vibration to 
total discomfort. For example, equations (20) and (22) would indicate to the ride- 
quality engineer whether a ride-quality problem derives from noise, from vibration, 
or from both. Additional information can be obtained from equation (21) to determine 
the discomfort contributed by individual noise octave bands and from equations ( 1 ) 
through (8) to determine discomfort contributed by individual vibration axes. Ulti- 
mately, the engineer could examine individual vibration and noise spectra to identify 
the particular spectral components that are the source of ride-quality problems and, 
then, institute remedial measures. Tht. potential of this approach in determining 
comfort trade-offs due to noise and vibration is apparent. Such a comprehensive 
ride-quality design and assessment tool has heretofore been unavailable. 
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DESCRIPTION OF aJMPUTER PROGRAM 
G«nf*r.il DHRorlptlon 

The ride comfort model dencrlhed In the i)recedinq sectlona Viaa i>een Implemented 
in a uner-orlented cnm|)iiter proqram, named RIDEyiJL, on the NASA Lanqley Research 
Center computer system. This proqram, which is available throuqh COSMIC (ref, 13), 
can be exercised by either of two basic approaches. One afiproach re<juires the user 
to analyze noise and vibration data to obtain the weiqhted rms acceleration levels 
and A-weiqhted octave-band noise levels required as model input. The resultant data 
are then input to the model in either a batch processlnq or an Interactive mode of 
operation. The second approach requires only analog tape recordings of the required 
data. These recordings are digitized, calibrated, and processed through programs 
available at Langley for time series analysis (see ref. 14) and noise analysis. 
Outputs of these programs are provided in a file format suitable for input to the 
computer program. These files are accessible for either batch or interactive pro- 
cessing. Thus, this approach requires no pre-analysis of the physical data by the 
user. 


The computer program RIDEQUL is written in Control Data Corporation's FORTRAN 
Version 5, based on FORTRAN 77, for the CDC CYBER 170 series computers operating 
under the Networlt Operating System at the NASA Langley Research Center, it utilizes 
one specific library subroutine available at Langley and the NAMELIST method for 
input and, therefore, may not be directly compatible with other systems even with the 
subroutine present. However, only minor adjustments should be required to run the 
program anywhere that a compiler for FORTRAN 77 is available. The data structures 
are those in use at Langley and may need to be adapted to the user's computer system. 


General Program Requirements 

To run the RIDEQUL program requires only that vibration, duration, and noise 
data be available for each situation to be analyzed. However, there are options for 
how these data may be available. There must be vibration present in at least one 
axis for the program to analyze a situation. This vibration may be characterized as 
random or sinusoidal, but not both. Random vibration is relevant in all five axes 
(vertical, lateral, longitudinal, roll, and pitch). Sinusoidal vibration is relevant 
in only three axes: vertical, lateral, and roll. The duration of the situation may 
or may not be ta)cen into account. Likewise, any noise present in the situation may 
or may not be taken into account. All input to the program should be nonnegative, as 
negative data are undefined. 

The vibration and the noise data are either input as specific physical param- 
eters known from prior analysis or input from TlFT-formatted (time history Interface 
file tape) multifile files obtained from the time series analysis or the noise analy- 
sis programs. The TIFT format for multifile files is described in appendix A. Each 
component file of a multifile file contains data for one situation. (These component 
files are Interchangeably referred to as files and as serials.) Data files in this 

format are obtained at the NASA Langley Research Center via the Central Dynamic Data 
Transcription Subsystem. 


R 


ORiGrr^AL i A;.:-: z; 

OF POOR (AlA..liY 

Proqriim InputH 

Vli)r .it~i on.- t!ih firnt input: nufidud In a apeci f ioatl on an t:o wh«th«r the vibra- 
tion 1 h random or nlnunoldal. If vibration In random, a TIKT-tormatted multifile 
fil»» containinq power apeotral denalty (jmd) levels for tlie axen of vibration in 
exfjected. All five axes (vertical, lateral, lonqitudlnal, roll, and pitr-h) may lx* 
Involved. The multifile file nhould lx? known locally an TAPK1 . Input In then the 
nerial numbera tor the tllea for each of the axea. Serial number 0 Indicates tl»at 
there la no vibration In an axis. The frequency and the pad level muat be chan- 
nela 1 and 2, rea|)ect. I vely, of the file, frequency taklnq the place of time In the 
format. The special case of a aituatlon with essentially purely sinusoidal vibra- 
tion la handled by indlcatlnq Oiat the vibration Is not random, vibration In only 
the vertical, lateral, and roll axes may iie specified to be sinusoidal and must then 
have been analyzed to determine the frequencies and root-mean-square (rma) accelera- 
tion levels associated with the dominant spectral components within each axis. Up 
to five components per "sinusoidal" axis may be used, although usually fewer will be 
required. The user chooses which spectral components to employ. 

Duration .- Duration input must be given in minutes but may be given as 0.0 if 
duration effect is not to be taken into account. 

Noise »- Whichever kind of vibration is present, noise is either input from a 
file or directly input as octave-band A-weighted sound pressure levels, if the user 
indicates tliat noise data are to come from a file, a TIPT-formatted multifile file 
known locally as TAPES is expected. The serial number of the appropriate file is 
then the noise input. This file must have the unweighted sound pressure levels for 
the one-third octaves beginning at 3.2 Hz as channels 2 through 31 of the data 
record. Extra one-third-octave sound pressure levels are ignored. Only one data 
record is expected in this file. If there is no noise to account for, a serial num- 
ber of 0 may be given, if it is indicated that the noise data are not on a file, 

then A-weighted sound pressure levels for the six relevant octave bands may be input 
directly. 


Program Execution 

The program may be executed in batch mode or from an interactive terminal. When 
accessed from the Langley computer complex directly, the compiled program may be 
found on an indirect access file named RQM under user number 036623C. The subroutine 
GETJO from ^e library UTLIB5, an indirect access file under the user number UTIL, 
must be available for loading at execution. The following commands will then invoke 
the program: 


FETCH , RQM/UN»0 366 2 3C . 
FETCH , UTLl B5 /UN^OTIL . 
LIBRARY, UTL1B5. 

RQM. 


If tlie program is obtained through aiSMlC (ref. 13), refer to the information 
that comes with the program and to the compilation and linkinci rctptirempnts of the 
local computer installation. 
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Mul ij)le a.turtHnnrt may te <turln<j tli« same i)ro.iram 

hH.jinnin.} of «x«out ion t he program roiaiimee a l l the ueer-def i nable 
either zero or true, depending on their type. 


exetniM ori 
Vfir i rthlen 


At the 
to he 


Uie pro.}ram Ik. run in the liatoh mode, a card deok may lx* read 
111 or a file of card ImaMeu may lx* nubmitted to the input queue. The input file con- 
uUtH of a Herien of NAMKMST input <!roui>M, each followed. If the vibration In slnu- 
Holdal by a card Imaqe of up to HO , haractern ui.ed to identify the Mltuation. If 
the vibration in random, the identification in taken from the header record on TAPR1 
for the firnt axin with vibration. Kach NAMKI.IST input .,roup provides the informa- 
tion needed for a jiartlcular situation. The NAMKI.IST conventions apply an<i are 
explained in Chapter *1 of the Cix: KOHTHAN Version 5 Reference Manual (ref. IS) under 

'NAMRhIKT Innut /Output." The NAMRLIST (jroup name is DSCPRM.S, The variables In It 
are 


RANDOMV 


ISER(5) 


Ixiqical variable set to 

TRUE If vibration Is random 

FALSE If vibration is sinusoidal 
Initialized to TRUE. 


Integer array of serial numbers of files on multifile file TAPE1 con 
talning random vibration power spectral densities for 5 axes; 
ISER(l) For vertical axis 
ISER(2) For lateral ai'is 
ISER(3) For longitudinal axis 
ISER(4) For roll axis 
ISER(5) For pitch axis 
Initialized to all 0. 


NVIBPKSO) 


Integer array of the number of vibration spectral components which 
will be input directly for 3 axes relevant to sinusoidal vibration. 
NVIBPKS(I) For vertical axis 
NVIBPKS(2) For lateral axis 
NVIBPKS(3) For roll axis 
Initialized to all 0. 


PRpPKS(5,3) Real array of frequencies, in hertz, of up to 5 vibration spectral 
components in each of 3 axes relevant to sinusoidal vibration. 
Initialized to all 0.0. 


RMSVBPK(5,3) Real array of rms levels, in g units or rad/sec^, of vibration spec 
tiv..l components corresponding to frequencies in FRQPKS. 
initialized to all 0.0. 


□RATION Real length of time, in minutes, of the situation. 

Initialized to 0.0, 


NOISAVD Logical variable set to 

TRUE If unweighted one-third-octave sound pressure levels are on 
multifile file TAPES 

FALSE If A-weighted octave sound pressure levels are to be directly 
entered at relevant center frequencies 
Initialized to TRUE. 
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Int^*<|er fu^rirtl ol Ihf* 1 i Ic on miiltitil<» til** TAiM'^H ooiii .» i ni n<t 

oriH-t hi ni-ool n VM Hoiiiui fironMiir<» 

I ni ti a 1 i 1o 0, 

H»*a I .»rray of fii»^ A-Wf*i (jtif »-*(i or-tavo fiontul |ir<vi;nir<* at tin* 

(t OHnter f r*tqut*fir! i »*s , 

Ini Llal i/atil to al l 0.0. 

Tin* NAMKl. I .S'l’ variahlftft .»rc not rni n i t i .t 1 i znO lx»twnr*n aitalysos of ni tiiat ' oi^, so 
all muHt be renpecified except those that are not different for the iu*w situation. 
However, if HANOOMV is true, NVIHOKS, KHpPK.S, and HMSVHI'K are iqnored and need tint lie 
reset. If KANDOMV is false, IHKR is iijtifired and need lutt ix* r<*set. If NOI.SAVn is 
true, ORNOISK need not he reset. If NOI.SAVD is false, NOISKH need nut )>e reset. See 
appendixes H and c for examples of hatch execution setups. 


NOI.Sint 


DHNOI.SKlh) 


Interactive mode .- If the proqram is being executed interactively, it asks for 
each piece of needed information before analyzing a situation. All numeric input is 
"list directed." This type of input is described in 'Chapter 5 of the CPC PORTRAN 
Version 5 Reference Manual, (ref. 15) under "List Directed Input." The A descriptc'r 
is used to '.nput character data so that keying in apostrophes around the characters 
is unnecessary. This descriptor is explained in the same chapter under "Format 
Specification," "edit Descriptors." Any question which contains "(Y/N)" requires the 
character Y to be entered for a yes response. Any other entry, including the charac- 
ter H or a carriage return without an entry, is taken to be a no rea-onut. Only 
initially does the program ask whether the vibration is random or :;'‘nusoidal and 
whether the noise sound pressure levels have been saved; these art; assumed to be the 
same for all situations during an interactive execution of the projraia. If the 
vibration is random, there are prompts to supply the serial numbers on multifile file 
TAPE1 for the files containing the vibration data. If the vibration is sinusoidal, 
there are prompts for supplying the frequency and rras acceleration levels for each 
sinusoidal component for each axis and for supplying the identification for the 
situation. If sound pressure levels are on multifile file TAPES, the corresponding 
serial number is requested. If not, the A-weighted octave sound pressure levels are 
requested. See appendixes B and C for examples of interactive program execution. 

(The vtser's input is seen as lower case.) After analyzing each situation, the pro- 
gram inquires, "Do you have more data to analyze?" to determine whether another 
situation is to be evaluated. 


Program Outputs 

Output from the program consists of a list of the frequency ranges over which 
the model is valid and, for each situation analyzed, the echoed input f)aram<'t» t s with 
a summary of the discomfort analysis. If the program is being executed inicr.ic- 
ti.vely, the list of ranges is optional. The summary of tlie discomfort analysis con- 
tains first the Identification for the situation, dt rived from either the vibration 
data file or from input. Then the discomfort components for each axis of vibration 
tf.re printed. If the vibration is random, the unweighted rms acceleration level of 
vibcatlon at all relevant frequencies, the weighted rms acceleration level of vibra- 
tion (>t these frequencies, and the discomfort level due to vibration in this axis are 
reported for each axis, if the vibration is sinusoidal, the discomfort levels due to 
each vibration spectral component and the discomfort level due to all vibration in 
each axis are reported. The combined-axes discomfort levels for the combination of 
the vertical, lateral, and roll axes and for the combination of the longitudinal and 
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con^^ibatiorto factor to l>e considered, the 

vihrlir a discomfo, n level due to trip duration and the duration-corrected 

vibration discomfort level are printed, if there ia any noise input, the contrtbu! 
tion to discomfort due to noise in eacli of the relevant octave bands is printed 
followed by Uie contribution to the <iiscomfort level due to noise in all relevant 
octaves and the noise-corrected, duration-correcte<i vibration discomfort. Finally 
the overall discomfort level is output. All values referred to as discomfort levels 
aie qiven in the subjective DISC units. See examples in the section "Application of 
Computer Proqram." All output is also on a local file named TAPE? so that additional 


SAMPLE APPLICATIONS 

section illustrates the application of the simplified NASA ride comfort 
model to ride environments measured on helicopters and automobiles. For each vehicle 
the various discomfort estimates are first obtained by manual computation using the 

paper and then by use of the computer program RIDEQUL. 

For these examples the vibration input is random. 

For purposes of illustration the reader should assume that time histories of 
vertical vibration and interior noise have been measured and recSSed duSng orufse 
aill^ T ^ currently operational helicopter. These measurements are available in 
magnetic tape. These tapes also contain sufficient calibration data 
lar^nalL^taoer^°" tape voltage output to appropriate engineering units. Simi- 
flocr of^ are available to define the vibration environment measured at the 

floor of an automobile. Three axes of vibration (vertical, lateral, roll) were 

Thriot^rdf^^ recorded in the automobile, but no interior noise data were collected. 
Ihe total discomf ore level associated with each vehicle is to be estimated as well 
as individual components of discomfort due to noise and vibration. 


Manual Computation 

Assume that the analog tapes me<> ured during flight of the helicop- 
ter w^re analyzed by the user to determine the parameters required as model input. 

For this case only vertical vibration and interior noise were measured. The 
unweighted rms acceleration, weighted rms acceleration, and A-weighted sound pressure 
levels within the six required octave bands are 


^^uwWert 0e1142g 
(^w^vejTt ” 0»0751g 
^ A ,63 ^^*8 

^A,125 ^ 81 *2 
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^A,2S0 = 


^'A,500 = ^'*•5 


*^A,1K = 88«1 

*^A,2K = 73.4 


where (^uw^vert measured unweighted vertical 
^A,i the A-weighted sound pressure levels in the 
center frequencies. The computation procedure is as 


rms acceleration 
octave bands with 
follows ! 


level and 
the indicated 


1. For vertical vibration equation (1) gives a single-axis discomfort of 


= 0.241 + 44.672(0.0751) 
= 3.60 DISC 


=•1 


The discomfort due to vibration in the other four axes is assumed to be zero. 

equation^ "(It measured, it is not necessary to apply 

to th, Lponlw t^« 


'"comb = °vert = 


desired without corrections for duration; 
and the total vibration discomfort from equation (20) is 


^^®dur “ 0 


Vb “ 

vibration discomfort is used in the computation of noise 

tion to t t^i follow. It also represents the contribution of vibra- 

tion to total estimated discomfort and is used again in equation (23). 

individufr^taw^teLs^orL^’’''^'®!’'” compute the noise discomfort produced by 

in.eger value prior to entering table V to obtain the coefficients of equation (21). 

a. For the 63-Hz octave band, !.» = 70, so Aa =; 0.7452 and R. - n 
(from table V) and 1 *470 (from table VI)* Thus, 3 


^N(63,70) ~ (0.7452 - 0.2337(3.60)1(1.470) 
= "0.14 ni.sc 
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since Di 


is negative, sat 


N(63,70) 
’’n(G3,70) 


b. For the 125-Hz octave band, = 81, so Aj = 1.9605, = -0.4704, 


Wi = 0.963. Thus, 


and 


^N(125,81) ~ M *9605 - 0.4704(3.60)1(0.963) 
= 0.26 DISC 


c. For the 250~Hz octave band, = 90, so Aj 
W i = 0.786. Thus, ^ 


3.2968, Bj = -0.6547, and 


^N(250,90) ” (3.2968 — 0.6547(3.60)1(0.786) 
» 0.74 DISC 


d. For the 500-Hz octave band, = 86, so A^, = 2.6649, - -0.5738, and 

Wi » 0.646, Thus, ■' J 


D 


N( 500,86) 


= (2.6649 - 0.5738(3.60)1(0.646) 
= 0,39 DISC 


e. For the 1000-Hz octave band, L« = 88, so Aj 
Wi = 0.688. Thus, ^ 


2.9732, Bj = -0,6145, and 


D 


N( 1000,88) 


= (2.9732 - 0.6145(3.60)1(0.688) 
= 0.52 DISC 


f. For the 2000-Hz octave band, L* = 73, so A- 
Wi = 1.448. Thus, ^ 


1.0312, Bj = -0.2995, and 


^N(2000,73) ” - 0.2995(3.60)1(1.448) 

= -0.07 DISC 

Since 2000,73) negative, set 
%( 2000,73) * ° 
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follows: ‘discomfort due to noise is computed by application of equation (22) aa 


N(i, j )max 


0.74 rise 


I o 


N(i,j) 


0 + 0.26 + 0.74 + 0.39 + 0.52 + 0 


= 1.91 DISC 


Thus, 


°Noise “ * 0.3(1.91 - 0.74) 

= 1.09 DISC 


6. ^tal estimated discomfort is computed by adding the vibration and noise 
components of discomfort using equation (23): 


s 1,09 + 3.60 
“ 4.69 DISC 


^tomobile .- Analog tapes of automobile vibration in the vertical lateral an<^ 

^ioi to Obtain tL wLgh?ed rL locale 

tion levels in each of these axes. Interior noise was not measured. 

in Se^piesencror the computational steps required to estimate passenger 'discomfort 

SLCcUratlo. levels “ai^htea ms 


Axis 

Unweighted rms 

Weighted rms 

acceleration 

acceleration 

Vertical 

0.0692g 

0.0405g 

Lateral 

0.0230g 

0.0153g 

Roll 

0.2454 rad/sec^ 

0.1943 rad/sec^ 


tion lre’’wua?foir^n computing vertical, lateral, and roll accelera- 

are equations (1), (2), and (4), respectively. Applying these equations gives 


= 0.241 + 44.672(0.0405) 
= 2.05 DISC 
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= 0.393 + 47.494(0.0153) 
“ 1.12 DISC 


i 

1 


i 

f 

I 


= -0.21 + 4.506(0.1943) 
= 0.67 DISC 


2. Since only vertical, 
sary to apply only equations 
follows : 


lateral, and roll vibrations were measured, it is neces- 
(9) through (13) to compute Dy^R. The procedure is as 


j a. Ran)c the single-axis discomfort components such that = 2 05 

[ D2 = 1.12, and D3 « 0.67. ’ ' 

b. Compute 

f 

J = 2.05/1.12 

t - 1.83 

t 

= \1(1 .12)^ + (0.67)^ 

» 1.305 

» 'J(2.05)^ + (1.12)^ + (0.67)^ 

= 2.430 


c. Since D^i > 0.88, then 


“combi = ♦ 1.65(2.43) 

« 3.57 

d. For this case, D4 > 0.40, so 
Vr ' “combi ■ ^-5’ 

e. Since longitudinal and pitch vibration were not measured, equation (18) 


D 


comb 


D.. 


VLR 


« 3.57 DISC 


and since vibration discomfort is not to be corrected for duration effect 
tlon (20) becomes ' 


equa- 


“viB 


cumb 


^3.57 DISC 




* 


i 

I 

I 

I 

I 

j 

i 

I 

) 


0 . 


■i. 
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vehicle interior noise was not of interest. Thus, in equa- 
KZi), Uuoise set to zero. Therefore, total estimated discomfort is given 


D 

T 


VI B 


3.57 DISC 


Application of Computer Program 

computer program RIDEQUL at NASA Langley Research Center the 
magnetic tape recordings that contain sufficient calibration 
Le SuL cau^^ engineering units. These magneUc t^es 

producf wra-£ormaSrf?i'e Processed through the time series analysis program to 

P oduce TIFT formatted files containing power spectral densities for the vibrations 

of interest and through the noise analysis program to produce similar files contain- 
ing unweighted noise levels in one-third octavL. Bot? sets of ?Ues a« «sS J^^he 
xecution of RIDEQUL. The A-weighting and summation of the one-third octaves into 
-weighted octave sound pressure levels is performed by RIDEQUL# 

namelist input, and sample output for batch 
processing of the helicopter example are presented in appendix B. The dialoque from 

UoU^Unet"" "r.^® i3 alsrpresen;ed! ThrS!g^^:d 

(solid line) and weighted (dashed line) vertical vibration spectra for the heliSptL 
are shown in figure 5, and the interior noise spectrum is given in fianlTe IZ lll 

Seiqhtino finctl because of the effect of the human comfort LnsiJvity 

weighting function for vertical vibration. cavity 

Automobile .- A sample run deck, the NAMELIST input, and sample output for the 
u ©mobile example, without one-third-octave noise data saved on a file, are pr«- 
sent®d in appendix C. The dialogue from an interactive session for the exampiria 

^ absence of noise may^be indi- 

therf "«">ber 0 or by typing an "n" to the question of whethf^ 

there is noise Present. The unweighted and weighted acceleration spectra for verti- 
cal, lateral, and roll vibration are shown in figure 7. 


CONCLUDING REMARKS 

from P®P®>=‘ Ja® presented and discussed a simplified ride comfort model derived 

DlifiLf?on developed by NASA Langley Research Center. The major sim- 

Plification involved frequency weighting of the vibration spectra of each ZIT in 

accordance with the respective human comfort sensitivity characteristics The mof- 

resulting spectrum provided the basic vibra-” 
uieh 1 ^ simplified model. Treatment of the noise in the moLl^was 

es3ma?S\o cLmS'uo^ procedures required to obtain discomfort 

bv examolee tev^n f ^ vibration environments were presented and illustrated 

tLn of^^lr automobile ride quality studies, implementa- 

tion of the simplified ride comfort model in a user-oriented computer program was 
diacuRs*.d, and applications of the program to helicopter and automobile ride environ- 
-nts w-rA also presented. This computer program, in conjunction with time series 
'Sis and noise analysis program capability, permits routine processing and con- 
m of physical ride environment data from analog magnetic tape recordings to 
fort estimates that relate directly to passenger acceptance of the environ- 
These discomfort estimates (model outputs) are in the form of single numbers 
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reflecting total estimated passenger discomfort, noise and vibration contributions to 
total estimated discomfort, discomfoi.t attributed to individual axes of vibration, 
and discomfort attributed to individual noise octave bands. 

The model presented herein has several important features that distinguish it 
from other methods of evaluating vehicle ride quality. First, the model has the 
unique capability of transforming individual elements of noise and vibration environ- 
ments into subjective discomfort units and then combining the subjective units to 
produce estimates of total discomfort as well as the various individual indices of 
discomfort. Secondly, the model is sensitive to changes in physical stimulus param- 
eters such as vibration frequency, vibration acceleration level, noise octave-band 
frequency, and noise level. This renders the model very useful for determining ride- 
quality design trade-offs and especially for assessing comparative vehicle ride com- 
fort and diagnosing the source of ride comfort problems. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 12, 1984 
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TIFT PILE FORMAT 

The TIFT (time history interface file tape) file format is based on the natural 
structure of time series data. Other types of data may also lend themselves to this 
structure. Data ate arranged into groups of each )iind of data, one group for each 
occurrence. Therefore, each "sample" (occurrence) represents the state of each 
"channel" (kind of recorded data) at a particular point in time. The time of the 
occurrence is the first channel and is counted when considering how many channels are 
recorded or how many words written. Because of the nature of time series, it is 
assumed that the values of the data in the first channel are in increasing order as 
one reads the file. 

This format structures each file with a header logical record followed by data 
logical records, each of the latter representing a sample (or an occurrence). The 
header record contains the integer serial number of the file, the integer number of 
channels of data in the file, the corresponding names of the channels of data in 
words of 10 characters each, the corresponding units of the channels of data in words 
of 10 characters each, and 8 words of 10 characters each that identify the data in 
the file. The data records contain one word for each channel of data. Each file is 
terminated with an end of file mark. A multifile file is terminated with an addi- 
tional end of file mark, resulting in two in a row. 

The TIPT format may be visualized as 
Record 1 - Header record 

Word 1 - Serial number (integer) 

Word 2 - Number of channels N (integer) 

Word 3 - Alphanumeric name of first channel, usually 'TIME ' 

Word 4 - Alphanumeric name of second channel 


Word N+2 - Alphanumeric name of Nth channel 

Word N+3 - Alphanumeric unit for first channel, usually 'SECOND ' 
Word N+4 - Alphanumeric unit for second channel 


Words 2N+2 - Alphanumeric unit for Nth channel 
Words 2N+3 to 2N+10 - Alphanumeric header identifying data 
rds 2 to n - Data records 

1 - Value for first channel, usually time value 

2 to N - Data values 
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II^PUT AND OUTPUT FOR BATCH OR INTERACTIVE PROCESSING OP 

HELICOPTER EXAMPLE 

Tlift followinq run deck and NAMELIST input result in hatch processinq by proqram 
RIDEQUL of the helicopter example qiven in the text. The output of this batch run is 
then presented. Finally an interactive session is presented in which the helicopter 
example is analyzed. 


Run Deck 


PTHFOl'l »T10C0»Cr*70000, 

MSFR»oroooox»yxxxxx)f. 

CHAPfiF,000000»l PC. 

OPLIVFP.XVX xxx rOUP PFLIVFPY 
FPTCM»TAPEl. 

FPTCH»TAPFP. 

FPTCH,pOM, 

FFTCM,l'TirR5/l'N«UTU. 

LI»PAPYU»TLIR5) 

MAPfOPP. 

POM, 

PFWINP,TAPEF»TAPE7, 

COPYrF»TAPF6»0UTPl*T. 

C0PYCFjiTAPP7»PUTPltT. 

FXIT. 


NAMELIST Input 



♦DSCPPHS I5FR»1*C#0»0»0» NOISER»0» 0RATinN«0.?5 $ 
POFCPP^S I?FP-2#C»0»0»C F 


P05CPPMS ISFP-3#0»C*C»0 t 
tnSCPPHF ISPP«A#C,0»C»0 * 
PHSCPPf^F ISFR«5»0,0»0»0 « 
♦nSCPPf^S IFFP«6»C#C#C»0 P 
tnscppi'F i<pp«7,o»o*o#c f 
♦Dscpp^s iSFp«?*o»c»o»r * 
»0SCPPMS I?FP*Q*('»C»C»0 t 


pnscpPRp iFPP*io#o»o#o»o t 
pnscpp*-s iSFP»n»o»o»o»o» npisfr-p 

fnSCPPMS ISFP«12»0,0»0.0» NOlSFR-3 
fPSCPRf^S ISFP«13»C»0#0»0» NPISFR-A 
♦Dscppr*^ i?FP*iA#o»o#o»o» KniSFP *5 
tOSrPRMC IFFR*15»0»0»C*0» NPI5FP«6 


% 

« 

% 

f 


■i 
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The followinq In part of the output resultinq from t)ie aht>ve run deck and 
NAMRhIST input; 


♦ ♦ OISCnMFOJ^T PROGRAh ♦♦ 16.19.3A. 01/aA/H4 

I 


IF Trtf VIBRATIIIN IS KANOOI# ALl FIVE AXFS ARE USEOf 


AXIS 

FREQUENCY 

RANGE 

VFRTICAl 

l.OO " 

30.00 

LATERAL 

1.00 " 

lO.CO 

longitudinal 

.50 - 

10.00 

ROLL 

.50 - 

5.50 

PITCH 

.50 - 

10.00 


IF THfc VIBRATION IS SINUSOIDAL# THF FOLLOWING THREE AXFS A#E USEOI j 

'I 

AXIS FREQUENCY RANGE 


VERTICAL 

LATERAL 

ROLL 


1.00 - 30.00 
1.00 - 10.00 
1.00 - A. 00 


CORRECTION FOR THE DURATION OF THE SITUATION CAN ACCOUNT FOR 
BETWEEN 1 AND 120 MINUTES. 

NOISE FOR OCTAVE BANOS WITH THE FOLLOWING CENTER FREQUENCIES WILL 
BE ACCOUNTED FOR IF THE NOISE LEVEL IS ABOVE 65 0B(A). NOISE ABOVE 
100 DB(AI CAN ONLY BE TREATED AS 100 DB(A). 


63. H2 
125. HZ 
250. HZ 
500. HZ 
1000. HZ 
2000. HZ 


SDSCPRMS 
RANOOMV • T» 

ISER * 11* 0# 0# 0# Of 

NVIBPKS • 0* 0# 0» 

FROPKS • 0.0# 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 

RMSVBPK • 0.0* O.D* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 0.0* 

0RA11UN • .25E+00* 

NOISAVO • T* 

NOISER “ ?* 

OBNOISE • 0.0* 0.0* 0.0* O.O* 0.0* 0.0* 

SEND 


* 


I 

I 

j 

1 

i 

I 

I 

1 

t 

I 

i 
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HtH SlM/NAVr-l/VIB/LEATHF*»WOOO 


TfST Nn. P-1 


RUN NT. 1 


FQH TMf; (/IPTICAl AXIS 

THE UNHl IGHTfD PIS IS .11A2 
THE PUGHTtU PMS IS .0751 
THl DISC'lMfURI out TO VIRPAIION 


IN THIS AXIS IS 3.60 DISC UNITS. 


COMBINED-AXIS OISCUMEDPT » OP VERTICAL* LATERAL* AND POLL AXES IS 3.60. 
COMBINED-AXIS DISCQNfORT FOP LONGITUDINAL AND PITCH AXES IS 0.00. 
TOE CQMltlNFO DISCOMFORT DUE TO VIBRATION IN ALL AXES IS 3.60. 


LENGTH OF TIMF LESS 

THAM 1 

MINUTt 

CANNOT BH ACCnUNTfc 

BE NU CURRLCTION FOR 

n^iF. 


DISCOMFORT FROM EACH 

OCTAVFl 


CENTER FREQUENCY 


nB(A) 

DISC UNITS 

63. 


69.8 

0.00 

129. 


8i«a 

.26 

290. 


89,8 

.74 

900. 



.39 

1000. 


88. 1 


2000. 


73.4 

0.00 

NOISE IN THE VIBRATION ENVIRONMENT 

CONTRIBUTES 1.09 1 


NOISE-CORRECTED* DURATION-CORRECTFO VIBRATION DISCOMFORT IS A. 69. 
THE OVERALL DISCOMFORT INDEX IS A. 69. 


Interactive Session 

The following is the dialogue resulting from interactive processing by program 
RIDEQUL. The user's responses to RIDEQUL's prompts are presented in lowercase 
letters. 


/ fetch, rqm, t apel , t ap©0 
FETCH COMPLETE. 

/fetch, ut 1 i b5/Mn=ut i I 
FETCH COMPLETE. 

/ 1 i brary ( ut 1 i bS) 
LIBRARY (UTLIB5) 

/map, off 
MAP, OFF. 

/rqm 


DISCOMFORT PRDSRAM ## 


15. 17. 34. 01/25/84 


THE PROGRAM YOU ftRE RUNNING WILL OflUCULftTE SUBJECTIVE DISCOMFORT 
INDICES FOR SITUATIONS FOR WHICH VIBRftTION AND NOISE LEVEL DATA ARE 
PROVIDED. WOULD YOU LIKE TO BEE THE LIMITS WITHIN WHICH THE 
CONTRIBUTIONS DUE TO VIBRATION AND NOISE LEVELS CAN BE ACCOUNTED FOR"? 
(Y/N) 

y 


IF THE VIBRATION IH RANDOM, ALl. FIVE AXES ARE USED: 


AXIS 


FREQUENCY RANGE 


VERTICAL 

LATERAL 

LONGITUDINAl, 

ROLL 

PITCH 


l .«0 ~ 30.00 
1.00 ~ 10. 00 
.50 10.00 

. 50 - 5. 50 

,50 - 10.00 
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0 


ir iHF viHRmiUN Hi niNuiiuiDni , thi rni i nwiNi. ihrm hxi s «Rr uiiM); 
nxin ^R^(,UJ^.N(:Y rhnui 


VERTICRU 
I ftlFRttL 
ROLL 


1. 00 10 . 00 
K 00 ^.00 


CURRECIIQN FOR THE UURH1 lUN UF IHh ‘illlJHIION CRN HCCOHNI M*R 
BETWEEN 1 HNU 1E*0 MINUTES. 

NPISF FOR OCTAVE BANDS WITH THE FOU UWINl-i CENTER ERrUUENClT li Wll- 
BE ACCOUNTED FOR IF THE NOIRE I EVEt IS ABOVE hU DB(A). NOJUf ABOVE 
100 DB(A) CAN QNIY BE TREATED AS 100 DB<A). 


'*'r 

03. 

HZ 


les. 

HZ 


£50. 

HZ 


500. 

HZ 


1000. 

HZ 


£000. 

HZ 


IS THE VIBRATION IN THESE SITUATIONS RANDOM (AS OPPOSED TO SINUSOIDAL)'^ 
(Y/N) 

? y 

IS THE NOISE DPTA SAVED ON A FILE (AS OPPOSED TO HAVINB TO BE KEYED IN 
AT THE TERMINAL)? (Y/N) 

? y 

GIVE THE SERIAL NUMBERS ON TAPEl FOR THE FIVE AXES’ ACflELEHATIONS ™ 

THE VERTICAL. LATERAL, LONGITUDINAL, ROLL, AND PITCH FORCES. (SERIAl 0 
WILL INDICATF NOT TO ACCOUNT FOR THAT PARTICULAR AXIS. AT LEAST ONF. 
SERIAL MUST BE NONZERO. ) 

? 1 1 , 0, 0, e, I2I 

WHAT IS THE LENGTH IN MINUTES OF THIS SITUATION? 

? 0.£5 

WHAT IS THE SERIAL NUMBER ON TAPES FOR THE NOISE DATA'' (0 WILL BE 
TAKEN TO MEAN NO NOISE CORRECTION IS DESIRED.) 

■> 2- 


SERIALS OF AXES’ ACCELERATIONS ARE 
VERTICAL 11 

LATERAL 0 

LONGITUDINAL 0 

ROLL 0 

PITCH 0 

TIME IN MINUTES IS . £S 
SERIAL FOR NOISE DATA IS £ 


HELI SlM/NAVY-l/VIB/LEATHERWUOD 


TEST NO. P-1 


RUN NO. 1 


FOR THE VERTICAL AXIS 
THE UNWEIGHTED RMS IR .114/: 

THE WEIGHTED RMB IS .0751 

THE DISCOMFORT DUE TO VIBRATION IN THIS AXIS IS ;L 60 DISC UNITS. 


COMBINED-AXIS DISCOMFORT FOR VERTK/At , LAU KOI , AND ROIL AXES IR S. f,0. 
COMBINED-AXIS DISCOMFORT FOR LPNRITUDINAI AND PITCH AXTB IR 0.00. 
THE COMBINED DISCOMFORT DUE TO VIBRATION IN Al l AXFR IS .-L(,0. 

LENGTH OF TIME LESS THAN 1 MINUTF CANNOT BF ACCOUNTFD FOR? IHl RF WII I 
BE NO CORRECTION FOR TIME. 




APPENDIX B 
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DlB^:0^rQRT f HUM hUUH UUIHVt * 


f FNItH rwlUUFNUY 

D&(D) 

G3. 

SS. 8 


81. 3 


8S. 8 

Mia. 

an, n 


88. 1 

3000. 

7:^. 4 

NniSt IN THF VIBHRTIQN 

ENVIRONMENT 

NOISE -CORRECTFD, DURATION CORRFCTE 


DIUC UNI IS 
1^. 

. 

. 7A 

• S» ‘ 

(A* 01^ 

CONTHlBUTEf; 1.109 DISCOMFORI UNITG 
I) VIHROTinN DIBCUMFORT I!i A. fc 9 . 


THE DVERHLL DISCOMFORT INDEX IS A.C9. 


DO YOU HOVE MORE DOTfl TO ONRLY/E (Y/N) 
"> n 


END OF PROGRAM EXECUTION. 

SUMMARY OF RESULTS IS ON TAPE?. 

0.805 CP SECONDS EXECUTION TIME, 
/route, tape?, dc-lp 
ROUTE COMPLETE. 



APPKNDIX I' 


1 ■ 


INHIT AND ODTHIT FOII MATCH OK INTKKACTIVK I'KOCKSS 1 NC OF 

AimiMOMIl.F. KXAMMl.F 

Thf* followlnq run rtiiH NAMJ'lif.'lT input; ruMult in Vi.iti-h pr<i<:*‘su i nq tiy |iroi|riim 

RIDKQIJD of tim aufomohi Iff f^xainjile qivcn in tho tnxt. Tlu* ouf|»it of tliiu hatcii run in 
th»*n |>r<‘!!<‘iil »‘<1 . Finally an i n( '■riuU ! v<* lU-UMion iu prfuu-nt oil irt wliii'li tlio antomoDiln 
example la analyKftd. 


Run Deck 


PTOFOUl ♦TlOOOfCr'^OOCO. 

i)SFR»coooooy»xxxxyxy. 
CMAPCfF» 0000 C 0 #LRC. 
nPlIVFP.VXV XXX vnilp rFLIVFRV 
FFTCH.TAPFI. 

FFTCH♦Pe^•. 

FFTCH»UTl Tn^/UNoUm , 
LTRPARVd'TLIPS) 

MAP, OFF. 

POM, 

PFWIMP,TAPF6,TAPF7. 

CnpvrP,TAPFf»OUTPfT. 

C 0 PVCF,TAPF 7 , 01 TPI^T. 

FXIT. 


NRHELIST Input 


^OSCPPMS ISFP*1»?»0#3#C# N'nispp«0» DPATTPN-O.O t 
♦OSCPPMS ISFP*A,«,C»F ,0 t 
tnSfPPMS I??Fpii7,P,0,P,P % 

«nsCPPMS lSFR« 10 ,n» 0 #l 7,0 f 
tDSCPP«S ISFP« 13 ,l<t, 0 , 15,0 « 
fnSCPPM^ IfFP«« 6 «, 5 ff 6 ,P, 5 Af ?,0 f 
UPSCPPM^ ISFR•»i^P,^ 6 P, 0 , 570 ,^ f 
♦nSCPRMfl ISFP- 57 ], 572 , 0 , 573,0 f 
tnSCPPMP ISFP- 57 A, 575 , 0 , 576,0 ^ 
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Output 


The followlnq is part of the output reaultinq from the above run deck and 
NAMKLTST input: 


♦ ♦ niSCUflFORT PROGRAM ♦♦ 09.3<5.<»0. 0l/?*5/RA 


IF THF VIHRATION IS KANOOMt AU FIVE AXES ARE USEOt 


AXIS 

FREQUfcNCY 

RANGE 

Vf RTICAt 

1.00 - 

30.00 

LATn«AU 

l.OO - 

10.00 

IQNGITUOINAl. 

.^0 - 

10.00 

ROLL 

.50 - 

5.50 

PITCH 

.50 - 

10.00 


IF THE VIBRATION IS SINUSOIDAL# THE FOLLOWING THREE AXES ARE USED* 
AXIS FREOUENCY RANGE 


VERTICAL 

LATERAL 

ROLL 


1»C0 30.00 

1. 00 - 10.00 

1.00 - 4.00 


CORRECTION FOR THE DURATION OF THE SITUATION CAN ACCOUNT FOR 
BETWEEN 1 AND 120 MINUTES. 


NOISE FOR octave BANOS WITH THE FOLLOWING CENTER PREOUENCIES WILL 
BE ACCOUNTED FOR IF THE NOISE LEVEL IS ABOVE 65 D8(A). NOISE ABOVE 
100 0B(A) -vMN ONLY BE TREATED AS 100 DBCA). 


63 . 

HZ 

125. 

HZ 

250. 

HZ 

500« 

HZ 

1000. 

HZ 

2000* 

HZ 


iOSCPRMS 






RANDOM V 

• 

T# 




ISER 

m 

566# 

569# 

Of 970f Of 

NVIBPKS 

m 

D# 0# 

0# 



PROPKS 

m 

0.0# 

0.0# 

0.0# 

0.0# 

RHSVOPK 

ifi 

0.0# 

0*0# 

o 

• 

o 

0.0# 

ORATION 

m 

0.0# 




NOISAVD 

m 

T# 




NQISPR 

« 

0# 




DBNOISF 

n 

0.0# 

0*0# 

0.0# 

u .0# 


0 . 0 # 

0 . 0 # 


0 . 0 # 


0 . 0 # 

0 . 0 # 


0 . 0 # 

0 . 0 # 


0 . 0 # 0 . 0 # 0 . 0 # 

O.C# 0.0# 0.0# 


0 . 0 # 

0 . 0 # 


0.0# 0.0# 0.0# 

OeO# 0.0# 0*0# 
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SEND 

NQ CORRECTION WILL BE MADF FUR TIMF. 
NO CORRECTION WILL BF MADE FOR NOISF 
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NASA/FORO TEST /LE ATHF RWO'TO 


test no, 16 MO, 10 


THE 

THE 

the 


FOR THE VERTICAL 
UNWEIGHTED RNS IS 
WEIGHTED RHS IS 
DISCQHFORT DUE TO 


AXIS 

.0692 

.0409 

VIBRATION 


IN THIS AXIS IS 2.09 DISC UNITS. 


FOR THE LATERAL 

The unweighted rns is 

THE WEIGHTED RMS IS 
THE DISCOMFORT DUE TO 

FOR THE ROLL 
the unweighted RMS IS 
THE WEIGHTED RMS IS 

the discomfort oue to 


AXIS 

.0230 

.0193 

VlflwATION IN THIS AXIS 

AXIS 

.2494 

.1443 

VIBRATION IN THIS AXIS 


IS 1.12 OISC UNITS. 


IS .67 DISC UNITS. 


THE COMBINED DISCOMFORT OUE TO VIORAtIon iHu Ix2s 


the OVERAlt DISCOMFORT INDEX IS 3.97, 


3.57. 

0 . 00 . 

3.57. 


Interactive Session 


/fetch, rqm, tapel 
FETCH COMPLETE, 
/fetch, ut 1 i b5/urt=ut i 1 
FETCH COMPLETE. 

/I ibrary (ut 1 ibS) 
LIBRPRY(UTLIB5) 

/map, off 
MPP, OFF. 

/rqm 


DISCOMFORT PROGRAM ** 


15.09. 21. 01/25/84 


CONTRIBUTIONS DUE TO VIBRRTION AND NOISE LEVELS CAN BE ACCOUNTED FOR? 


IF THE VIBRATION IS RANDOM, 
AXIS 

VERTICAL 
LATERAL 
LONGITUDINAL 
ROLL 
PITCH 


ALL FIVE AXES ARF USED: 

FREUUENCY RANGE 

1.CH9 - 3(0.010 
1.00 - 10,00 
.50 ~ 10.00 
.50 5,50 

.50 - 10.00 
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IF THE VIBRATION IS SINUSOIDAL, THE FOLLOWING THREE AXES ARE USEOs 


AXIS 


FREQUENCY RANGE 


VERTICAL 

lateral 

ROLL 


1.013 - 30.00 
1.00 10.00 
1. 00 - 4. 00 


CORRECTION FOR THE DURATION OF THE SITUATION CAN ACCOUNT FOR 
BETWEEN 1 AND 1£0 MINUTES. 


NOISE FOR OCTAVE BANDS WITH THE FOLLOWING CENTER FREQUENCIES WILL 
BE ACCOUNTED FOR IF THE NOISE LEVEL IS ABOVE 65 DB(A). NOISE ABOVE 
100 DB(A) CAN ONLY BE TREATED AS 100 Db<A). 


63. HZ 
laS. HZ 
£50. HZ 
500. HZ 
1000. HZ 
£000. HZ 


IS THE VIBRATION IN THESE SITUATIONS RANDOM (AS OPPOSED TO SINUSOIDAL)? 
(Y/N) 

’ y 

IS THE NOISE DATA SAVED ON A FILE (AS OPPOSED TO HAVING TO BE KEYED IN 
AT THE TERMINAL)? (Y/N) 

? n 


GIVE THE SERIAL NUMBERS ON TAPEl FOR THE FIVE AXES’ ACCELERATIONS -- 
THE VERTICAL, LATERAL, LONGITUDINAL, ROLL, AND PITCH FORCES. (SERIAL 0 
WILL INDICATE NOT TO ACCOUNT FOR TfiAT PARTICULAR AXIS. AT LEAST OWE 
SERIAL MUST BE NONZERO. ) 

? 568,569,0,570,0 

WHAT IS THE LENGTH IN MINUTES OF THIS SITUATION? 

? 0.0 

NO CORRECTION WILL BE MADE FOR TIME. 

IS THERE NOISE PRESENT? (Y/N) 

? n 


SERIALS OF AXES’ ACCELERATIONS ARE 
VERTICAL 56B 

LATERAL 569 

LONGITUDINAL 0 

ROLL 570 

PITCH 0 

TIME IN MINUTES IS 0.00 
NO CORRECTION WILL BE MADE FOR NOISE. 


NASA/FORD TEST MV -1 /LEATHERWOOD 


TEST NO. 16 


RUN NO. 10 


FOR THE VERTICAL AXIS 

THE UNWEIGHTED RMS IB . 069£ 

THE WEIGHTED RMS IB . 040^, 

THE DISCOMFORT DUE TO VIBRATION IN THIS AXIS IS £.05 DISC UNITS. 

FOR ThT LATERAL AXIS 

THE UNWEIGHTED RMB IS . 0£30 
THE WEIGHTED RMS IS .0153 

THE DISCOMFORT DUE TO VIBRATION IN THIS AXIS IS U l£ DISC t?NITS. 


FOR THE ROLL AXIS 

THE UNWEIGHTED RMB IS .£454 
THE WEIGHTED RMS IS ,1^4 -; 

THE DISCOMFORT DUF TO VIBRATION IN THIS AXIS IS 


.67 DISC UNITS. 
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COMBINED-AXIS DISCOMFORT f UR VERTICAL, LATERAL, AND ROl L AXE:S IS 
COMBINED-AXIS DISCOMFORT FOR LONGITUDINAL AND PITCH A>1 S IS 
THE COMBINED DISCOMFORT DUE TO VIBRATION IN ALL AXES IS 

THE OVERALL DISCOMFORT INDEX IS 3,57. 


DO YOU HAVE MORE DATA TO ANALYZE? (Y/N) 
? n 


END OF PROGRAM EXECUTION. 

SUMMARY OF RESULTS IS ON TAPE?. 

4.354 CP SECONDS EXECUTION TIME. 
/rout», tape?, dc=lp 
ROUTE COMPLETE. 


S u 
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TABLE I.- HUMAN COMPORT SENSITIVITY WEIGHTING FACTORS FOR EACH AXIS 

OP VIBRATION 


WEIGHTING FACTORS 


FREQUENCY# 


H7 

VERTICAL 

LATERAL 

RQU 

LONGITUDINAL 

PITCH 

.50 

.000 

*000 

.500 

.010 

• 60 

.000 

.000 

• 518 

• 030 

.70 

.000 

*000 

.536 

.070 

• 80 

.000 

*000 

.556 

.170 

.PO 

*000 

*0C0 

.572 

• 610 

1.00 

• 610 

.600 

• 590 

• 890 

1.25 

• 618 

.700 

.693 

1^000 

?.00 

• 660 

1.000 

1.000 

1^000 

3 .00 

• 660 

*750 

*670 

l^OOO 

^.00 

• 830 

.600 

• 660 

1^000 

<».75 

.890 

.525 

• 653 

UOOO 

5.00 

«910 

,500 

.650 

• 890 

5.50 

.^55 

• 665 

• 660 

• 650 

6.00 

1.000 

*630 

• 000 

• 230 

O.50 

• 925 

*600 

.000 

• 130 

f.OO 

.850 

.370 

.000 

• 080 

7.50 

.705 

.350 

.000 

• 050 

S.OC 

.760 

.33C 

,000 

• 060 

9.00 

.710 

*300 

*000 

• 020 

10.00 

.660 

*250 

.000 

• 010 

11.00 

.660 

*C00 

.000 

• 000 

12.00 

.620 

*000 

• 000 

.000 

n.oo 

.580 

.000 

.000 

• 000 

lA.OO 

.5^0 

*000 

.000 

• 000 

15.00 

.500 

*000 

.000 

.000 

16.00 

*670 

*000 

.000 

• 000 

17.00 

*660 

*000 

• 00 c 

• 000 

16.00 

*600 

*eoo 

*000 

.000 

19,00 

.380 

.000 

.000 

.000 

?0.00 

.370 

*000 

.000 

• 000 

?1.00 

.350 

.000 

• 00 c 

• 000 

?2.00 

.360 

*000 

*000 

• 000 

23.00 

• 330 

*000 

*000 

• 000 

26.00 

.330 

• 000 

• 000 

.000 

25.00 

.300 

.000 

• 000 

• OtC 

26.00 

• 300 

.003 

.000 

• 000 

27.00 

• 300 

.000 

• 000 

• 000 

28.00 

.300 

.000 

• 000 

• 000 

29.00 

• 300 

*000 

,000 

• 000 

30,00 

*300 

rOOO 

.000 

• 000 
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TABLE V.- SLOPE AND INTERCEPT OOEPPICIENTS EOR EQUATION (21) 


Noise level. 
La 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 


Intercept, 
Aj, DISC 

Slope, 

Bj 

Noise level, 
La 

Intercept, 
A*j , DISC 

Slope, 

0.3447 

-0.1219 

83 

• 

2.2294 

-0.5118 

.4172 

-.1445 

B4 

2.3718 

-.5329 

.4935 

-.1669 

85 

2.5164 

-•5533 

.5736 

-.1893 

86 

2.6649 

-.5738 

.6575 

-•2116 

87 

2.8172 

-.5942 

.7452 

-.2337 

88 

2.9732 

-.6145 

.8368 

-.2558 

89 

3.1330 

-.6346 

.9320 

-.2777 

90 

3.2968 

-.6547 

1 .0312 

-.2995 

91 

3.4642 

-.6746 

1.1340 

-.3212 

92 

3.6354 

-.6944 

1 .2408 

-.3429 

93 

3.8104 

-.7142 

1.3512 

-.3644 

94 

3.9893 

-.7338 

1 .4654 

-.3858 

95 

4.1720 

-.7533 

1.5835 

-.4071 

96 

4.3574 

-.7724 

1 .7055 

-.4284 

97 

4.5486 

-.7921 

1.8311 

-.4494 

98 

4.7426 

-.8113 

1 .9605 

-.4704 

99 

4.9404 

-.8304 

2.0938 

-.4913 

1 

100 

5.1421 

-.8494 


TABLE VI OCTAVE-BAND WEIGHTING FACTORS 


Octave-band 

Weighting 

center 

factor, 

frequency, Hz 


63 

1 .470 

125 

.963 

250 

.786 

500 

.646 

1000 

.688 

2000 

1 .. . 

1 .448 
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(b) Lateral axis. 

Figure 1 Human comfort sensitivity weighting functions 

of vibration, 


for various axes 



FREQUENCY. HERTZ 


(c) Roll axis. 
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(d) Longitudinal and pitch axcss. 


Figure 1 Concluded 



Power spectral density gvHz 



Frequency. Hz 


Figure 2.- l^eal spectral characteristics used in the development 
of the simplified NASA ride comfort model. 



ORIGiNAL PA”?: M 
OF POOR QUALITY 


H 'I-' - 
( ! 

. ; - 
aj 
c J 


X 





0 I i • 1 L- 

0 .01 .0.? .03 .0<4 


o 

Q> 

O 



O 

O 

O 


O 


* * ‘ ‘ * t 

.05 .06 .07 .08 .09 .10 


Weighted rnis acceleration, g unit 


Figure 3,~ Predicted discomfort (ref» 12) vs. 
weighted rms vertical acceleration. 
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Figure 4.- Predicted discomfort (ref. 12) vs, 
weighted rms lateral acceleration. 
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VIBRATION PSD, gVHz 
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Figure 5,- Unweighted and weighted vertical vibration spectra for 

the helicopter example. 




erior noise spectrum shape used in helicopter example 
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(b) Lateral axis. 
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(c) Roll axis. 


Figure 7,- Concluded 



